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Introduction 


*  Fuel  film  cooling  is  critical  for  high  performing  boost  engines 
using  the  Oxygen  Rich  Staged  Combustion  Cycle 

—  Reduces  heat  load  and  closes  the  power  balance 
—  Reduces  the  required  coolant  channel  pressure  drop 
—  Must  be  efficient  or  I5P  will  be  affected 
—  Film  cooling  schemes  are  not  optimized 
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Introduction 


*  FFC  inlets  can  have  complex  instabilities 


—  Acoustic  resonance  modes 
—  Kelvin-Helmholtz  instabilities 
—  Shear  layer  impingement 

*  Difficult  to  capture  computationally 

—  High  density  gradients 
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Carbon  Deposits 


—  Thermal  radiation,  soot  formation,  etc. 
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Previous  work 


*  Comparison  with  experiments  in  a 
subscale  fuel  film  cooled  thrust 
chamber 

—  2D  steady  state  RANS  with  ideal 
gas  and  equilibrium  chemistry 
(Himansu  et  al.  2014) 

—  Over-predicted  film  cooling 
effectiveness  near  FFC  inlet 

—  Under-predicted  performance 
further  downstream 

—  Does  not  capture  trends  in 
changing  mass  flow  rate 


Normalized  Distance  Downstream  of  FFC  Slot 


Hypothesis;  Unsteady,  pulsing,  flow  results  in  a  measured  heat  flux  that  is  the  average 
of  an  uncooled  and  an  overcooled  wall.  Steady  RANS  cannot  capture  this  effect. 
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Current  Objectives 


*  Capture  unsteady  flow  effects  in  FFC  simulations 

—  Determine  most  relevant  physics  to  thermal  management 
simulations 

—  Incorporate  complex  physics  into  simulation  suite 
—  Reproduce  results  seen  experimentally  in  wall  heat  flux 

*  Determine  if  added  cost  of  unsteady  simulations  improves 
predictive  capabilities 
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Numerical  Methods 


*  Using  General  Equation  and  Mesh  Solver  (GEMS) 

—  Unstructured  mesh  solver 

—  Second-order  accurate  in  time  and  space 

*  Ideal  gas  law  used  for  initial  study 

—  More  stable  calculations,  but  not  all  fluid  properties  are 
represented  accurately. 

—  Density  ratios  controlled  by  modifying  molecular  weight 
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Numerical  Methods 


*  Unsteady  time-accurate  simulations 

—  Detached  eddy  simulations  (DES) 

•  LES  subgrid  models  used  in  region  where  grid  is  fine 

•  Improved  performance  over  RANS  without  cost  of  full  LES 

•  Turbulence  computed  using  the  k-o)  model 

*  Steady  state  simulations 

—  Reynolds-Averaged  Navier-Stokes  (RANS) 

—  Used  as  initial  condition  for  unsteady  simulations 
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Baseline  and  Parametric  Cases 


*  All  2D,  same  baseline  geometry  as  Himansu  et  al. 

*  No  conjugate  heat  transfer 

—  Adiabatic  wall  temperature  calculations  in  place  of  wall 
heat  flux 


d 
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Parametric  Studies 


*  Baseline  case  chosen  from  in-house  experimental  setup 

*  Parametric  studies  will  vary  around  this  set  of  conditions 

*  Results  can  be  compared  to  experimental  trends 


Table  1.  Flow  conditions  for  baseline  case 


Pressure  (MPa) 

5.2 

Slot  Width  (m) 

5.6  X  10-'* 

FFC  Temperature  (K) 

350 

Main  chamber  temperature 

3680 

FFC  mass  flow  rate  (kg/s) 

0.41 

Main  chamber  mass  flow  rate 

25.64 

FFC  velocity  (m/s) 

1.1 

Main  chamber  velocity 

237.5 

FFC  density  (kg/m^) 

638 

Main  Chamber  density 

4.4 

FFC  Reynolds  number 

3.8  X  10'‘ 

Main  chamber  Reynolds  number 

2.4  X  10® 

Mass  flow  ratio 

1.6  X  10"2 

Mach  number 

0.2 

Momentum  ratio 

7.9  X  10”® 

Density  Ratio 

145 
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Parametric  Studies 


*  First  study:  Vary  main  chamber  mass  flow  rate,  FFC  mass 
flow  rate  to  produce  range  of  mass  flow  ratios  at  constant 
fuel  density 

—  Initial  study  to  determine  if  unsteadiness  is  significant 

*  Second  study:  Vary  density  of  fuel  film,  hold  momentum 
ratio  and  mass  flow  ratio  constant 

—  Determine  if  observed  frequencies  depend  on  fuel  density 

*  Third  Study:  Vary  slot  width,  hold  mass  flow  ratio, 
momentum  ratio,  and  density  constant 

—  Determine  if  velocity  trends  can  be  collapsed 
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Parametric  Studies 


Parametric  study 

Constant  density  and  slot 
width,  varying  mass  flow 
rates  and  momentum 

Constant  mass  flow  rates 
and  momentum,  varying 
density  and  slot  width 

Constant  density,  varying 
slot  width  and  momentum 

Reynolds  number 

«  10^  -  10® 

10® 

«  10® 

Slot  width  (m) 

5.6  X  10-4 

3.8-11x10-4 

lISJiWIiilkllibttHIII 

Main  chamber  density  (kg 
/m^) 

4.31  -  4.44 

4.42 

4.42 

Main  chamber  velocity 
(m/s) 

ii50m273NN 

233 

233 

Fuel  density  (kg  /  ni^) 

638 

146  -  432 

295 

Fuel  velocity  (ni/s) 

0.7  -  5.4 

1.2 -1.6 

lii^lillll 

Momentum  thickness  (m) 

1.5 -2.4  X  10-4 

1.6  X  10-4 

1.6  X  10-4 

Mass  flow  ratio 

0.01  -  0.075 

0.014 

0.014 

Density  ratio 

145 

33-97 

67 

Momentum  ratio 

10-4  _  ;L0-3 

«  10-4 

10-4 

Mach  number 

0.05-0.33 

0.2 

0.2 

Quantities  Varied 
in  Parametric  Study 
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Parametric  Studies 


Parametric  Study 

#1 

#2 

#3 

Reynolds  number 

0(10^  -  10«) 

O(IO^) 

O(IO^) 

Slot  width  (m) 

5.6  X  10“'* 

3.8-11x10-4 

3.8-11x10-4 

Main  chamber  density  (kg/m^) 

4.4 

4.4 

4.4 

Main  chamber  velocity  (m/s) 

50—273 

233 

233 

Fuel  density  (kg/m^) 

638 

146—432 

295 

Fuel  velocity  (m/s) 

0.7— 5.4 

«  1.4 

1—3 

Mass  flow  ratio  (kg/s) 

0.01—0.075 

0.014 

0.014 

Density  ratio 

145 

33—97 

67 

Momentum  ratio 

0(10“'*  -  10“^) 

0(10-4) 

0(10-4) 

Mach  number 

0.05—0.32 

0.20 

0.20 

.  ..  . 
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nt 
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Unsteady  DES 


*  Baseline  case 
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Comparison  of  PANS  and  DES 


*  Baseline  case 


Fuel  mass  fraction  Temperature 


Steady  RANS  I  Steady  RANS 
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T  emperature  ( K) 


Adiabatic  Wall  Temperatures 


*  Adiabatic  wall  temperature  profiles 


- ► 

Increasing  mass  flow  ratio 
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Adiabatic  Effectiveness 


*  Non-dimensionalize  as  adiabatic  effectiveness: 

T'f  Taw 

ri  =  7p - — 

‘f  ~  ^FFC 

Tf  =  freestream  flow  temperature 
=  adiabatic  wall  temperature 
Tpp(-  =  fuel  film  injection  temperature 

*  Allows  for  easy  comparison  of  all  three  parametric  studies 

*  Compare  to  experimental  trends  for  qualitative  analysis 
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Adiabatic  Effectiveness 


Mass  flow  varies,  constant  density/slot  width 


FFC  Velocity  - 

(m/s) 

Main 

Velocity 

50—371 

Main  Mach 

0.05—0.32 

638.51 

(kg/mO 

Main  Density 

4.41 

Density 

Ratio 

145 

5  6  X  lO--* 

width  (m) 

Mass  flow 
ratio 

1.0  X  10-2 
-  7.5  X  10-2 

Momentum 

ratio 

3.1  X  10“^ 

-  1.7  X  10-2 
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Adiabatic  Effectiveness 


Both  mass  flow  rates  vary,  constant  density/slot  width/mass  flow  ratio 

FFC  Velocity 
(m/s) 

0.7— 1.4 

Main  Velocity 

150—300 

Main  Mach 

0.13—0.26 

FFC  Density 
(kg/m^) 

640 

Main  Density 

4.4 

Density  Ratio 

145 

FFC  slot  width 
(m) 

5.6  X  10"^ 

Mass  flow  ratio 

1.6  X  10“^ 

Momentum 

ratio 

7.9  X  10"= 
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Adiabatic  Effectiveness 


0.4 


0.2 


0  H - 1 - 1 - 1 - 1 - 1 - 1 - 1 

0  0.01  0.02  0.03  0.04  0.05  0.06  0.07 

x(m) 


Density  varies,  slot  width  varies,  constant  mass  flow/momentum 


FFC  Velocity  .  . . 

(m/s) 

Main 

Velocity 

233.4 

Main  Mach  0.2 

FFC  Density 
(kg/mO 

Main  Density 

4.411 

Density  33_gg 

Ratio 

FFC  slot  3.8  X  10"'‘ 

Mass  flow 

1.4  X  10-2 

Momentum  ^2  x  lO""^ 
ratio 

width  (m)  —  1.1  X  10“^ 

ratio 
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Adiabatic  Effectiveness 


Slot  width  varies,  constant  density/mass  flow 


FFC  Velocity 
(m/s) 


Main 

Velocity 


233 


FFC  Density  _ 

(kg/m^) 

FFC  slot  3.8  X  10“^ 
width  (m)  —  1.1  X  10“^ 


Main  Density  4.4 

Mass  flow  4  X  10~^ 
ratio 


Main  Mach  0.2 

Density 

Ratio 

Momentum  6.1  x  10“^ 
ratio  —  1.8  x  10“"*^ 
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steady  v.  Unsteady  Effectiveness 


•  Main  mass  flow  rate  changes,  density  and  slot  width  constant 
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Adiabatic  Effectiveness 
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Power  Spectral  Density 


*  Baseline  case 

20  - 1 - 1 - r 


N 


0  5000  10000  15000  20000 


Frequency  (Hz) 
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Power  Spectral  Density 


*  First  parametric  study 

—  Increasing  main  chamber  mass  flow,  FFC  mass  flow  rate 
held  constant 


Frequency  (Hz)  Frequency  (Hz)  Frequency  (Hz) 


Increasing  mass  flow  ratio 
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Nondimensionalization 


*  Frequency  of  fuel  film  oscillations  is  not  constant 

—  Seems  to  switch  between  two  modes 


Helmholtz  number  for  each  case: 


He  = 


2nf 

cd 


f  =frequency 
c  =  speed  of  sound  (m/s) 
d  =  depth  of  FFC  inlet 

““  Integer  values  mean  inlet  is  functioning  as  resonator 
—  Calculated  He  for  three  most  dominant  frequencies  in  each  case 
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Helmholtz  number 


Helmholtz  number  v.  Mach  number 


*  Resonant  frequencies  for  first  parametric  study 


3.5 


A  ^  ^  A.aA  i 


2.5 


□ 


□ 


□  □ 


□ 


□  ■  □  □ 


□ 


□ 
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o 


o 
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♦  ♦  ♦  ♦  o 


o 


o 


0.5 


O  First  Resonance 
□  Second  Resonance 
AThird  Resonance 


0  0.05  0.1 


0.15  0.2 

Mach 


0.25  0.3  0.35 
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Conclusions 


*  Significant  differences  between  steady,  unsteady  simulations 

—  Capture  of  additional  flow  features  justifies  the  increased 
computational  cost  of  unsteady  simulations 

—  Unsteady  simulations  predict  lower  average  heat  flux  past  x/d=15, 
consistent  with  original  hypothesis 

*  Mass  flow  ratio  is  primary  driver  of  adiabatic  effectiveness 

—  Density  ratio  is  also  a  factor 

*  Unsteadiness  in  these  cases  driven  by  Helmholtz  mode  of  inlet  slot 

—  Inlet  geometry  important  for  film  cooling  effectiveness 

*  Still  significant  discrepancies  between  CFD  and  experimental  results 

—  May  be  due  to  soot  deposition,  radiative  heat  transfer,  or  other 
complex  physical  phenomena 
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Future  Work 


*  Three-dimensional  simulations 

—  Transverse  waves  may  affect  wall  temperature  profile 
—  May  also  model  more  realistic  FFC  inlets 

*  Real  Gas  Equation  of  State 

—  Capture  supercritical  properties 

*  Introduce  additional  physics 

—  Soot  formation  and  deposition 
—  Radiative  heat  transfer 
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